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ABSTRACT 

Two  dimensional  solitary  and  periodic  waves  in  water  of  finite  depth  are 

considered.  The  waves  propagate  under  the  combined  influence  of  gravity  and 

surface  tension.  The  flow,  the  surface  profile,  and  the  phase  velocity  are 

functions  of  the  amplitude  of  the  wave  and  the  parameters  l  »  A/H  and 
2 

t  -  T/pgH  .  Here  X  is  the  wavelength,  H  the  depth,  T  the  surface 
tension,  p  the  density  and  g  the  gravity.  For  large  values  of  1  and 
small  values  of  the  amplitude,  the  profile  of  the  wave  satisfies  the  Korteweg 
de  vries  equation  approximately.  However,  for  t  close  to  1/3  this  equation 
becomes  invalid,  in  the  present  paper  a  new  equation  valid  for  T  close  to 
1/3  is  obtained.  Moreover,  a  numerical  scheme  based  on  an  integro-differential 
equation  formulation  is  derived  to  solve  the  problem  in  the  fully  nonlinear 
case.  Accurate  solutions  for  periodic  and  solitary  waves  are  presented,  in 
addition,  the  limiting  configuration  for  large  amplitude  solitary  waves  when 
t  >  ^  is  found  analytically.  Graphs  of  the  results  are  included. 
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SIGNIFICANCE  AND  EXPLANATION 


The  study  of  the  influence  of  surface  tension  on  nonlinear  water  waves 
has  attracted  much  attention  in  recent  years  because  of  the  challenging 
mathematical  difficulties  associated  with  the  non-uniqueness  of  the 
solutions.  Several  numerical  schemes  are  now  available  to  compute  capillary 
gravity  waves  in  water  of  infinite  or  moderate  depth.  (Schwarts  and  Vanden- 
Broeck  (1979)/  Chen  and  Saffman  (1980)). 

In  the  present  work  we  present  a  numerical  scheme  which  enables  us  to 
compute  waves  in  the  limit  as  the  ratio  of  the  depth  versus  the  wavelength 
tends  to  zero.  A  problem  of  non-uniqueness  is  discovered  and  discussed. 
Moreover/  a  new  equation  analogous  to  the  KOrteweg  de  Vries  equation  is 
obtained.  In  addition,  the  limiting  configuration  of  steep  depression 
solitary  waves  is  obtained  analytically. 
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SOLITARY  AMO  PtRIODXC  GRAVITY-CAPILLARY  HAVES  OF  FINITE  AMPLITUDE 


J.  X.  Hunter  and  J.-M.  Vanden-Broeck 


1.  Introduction 

Approximate  solutions  for  gravity  solitary  and  cnoidal  waves  of  small  amplitude  ware 
obtained  by  Rayleigh  (1876),  and  Korteweg  and  de  Vries  (1895).  These  results  were  derived 
systematically  by  Keller  (1948)  who  calculated  a  first  order  perturbation  solution  in 
powers  of  the  wave  amplitude.  This  work  was  extended  to  second  order  by  Lei tone  (1960). 

More  recently,  accurate  fully  nonlinear  solutions  for  gravity  solitary  waves  were 
obtained  by  Longue t-Higgins  and  Fenton  (1974),  By ett- Smith  and  Ionguet-Higgina  (1976), 
Hitting  (1975)  and  Hunter  and  Vanden-Broeck  (1982).  A  review  of  some  of  these  theories  can 
be  found  in  Miles  (1980). 

Accurate  solutions  for  periodic  gravity  waves  in  water  of  finite  depth  were  obtained 
by  Schwarts  (1974),  Ookelet  (1977),  Vanden-Broeck  and  Schwarts  (1979),  and  Rienecker  and 
Fenton  (1981).  The  effect  of  surface  tension  on  periodic  waves  was  investigated  by  Crapper 
(1957),  Harrison  (1909),  Hilton  (1915),  Pearson  and  Fife  (1961),  Schwarts  and  Vanden-Broeck 
(1979)  and  Chen  and  Saffman  (1980).  For  a  review  of  these  calculations  see  Schwarts  end 
Fenton  (1982). 

The  effect  of  surface  tension  on  solitary  waves  was  first  considered  by  Korteweg  and 
de  Vries  ( 1895 ) .  They  discovered  that  depression  solitary  waves  can  exist  for  sufficiently 
large  values  of  the  surface  tension.  A  systematic  perturbation  calculation  was  attesipted 
by  Shinbrot  (1981).  However,  his  results  are  partially  incorrect  because  he  excluded  the 
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possibility  of  dsprsssion  mtis.  a  first  order  psrturbstion  solution  allowing  depression 
waves  was  derived  by  Vanden-Broeck  and  Shen  (1982)  and  Benjaain  (1982). 

This  perturbation  calculation  is  invalid  when 

t  -  ~~~~2  (1.1) 

pgH 

is  close  to  1/3.  In  (1.1)  T  is  the  surface  tension,  p  the  fluid  density,  g  the 
gravitational  acceleration,  and  H  is  the  undisturbed  depth  of  the  fluid. 

In  this  paper  we  present  a  perturbation  calculation  valid  near  t  -  1/3.  In  addition, 
we  solve  the  exact  nonlinear  equations  numerically. 

In  Section  2  we  formulate  the  nonlinear  problem  and  briefly  review  same  of  the 
classical  perturbation  calculations.  In  Sections  3  and  4  we  describe  the  perturbation 
calculation  valid  for  1  close  to  1/3.  In  Section  S  we  reformulate  the  problem  as  an 
integro-differential  equation  on  the  free  surface,  which  allows  us  to  calculate  solitary 
and  periodic  waves  of  arbitrary  amplitude.  A  numerical  scheme  for  solving  the  integro- 
differential  equation  is  presented  in  Section  6.  The  numerical  method  is  similar  in 
philosophy  if  not  in  details  to  the  scheme  derived  by  Vanden-Broeck  and  Schwartz  (1979). 

The  results  of  the  numerical  computation  are  discussed  in  Section  7.  In  addition  the 
limiting  configuration  for  large  amplitude  solitary  waves  when  t  >  1/2  is  found 
analytically  in  Section  7. 


2.  formulation  and  Classical  Perturbation  solution* 


W«  consider  a  two  dimensional  progressive  wave  in  an  irrotational  incompressible 
inviscid  fluid  having  a  free  surface  with  surface  tension  acting  upon  it,  and  bounded  below 
by  a  flat  horizontal  bottoa.  We  take  a  frame  of  reference  in  which  the  flow  is  steady, 
with  the  X-axis  parallel  to  the  bottoa,  and  with  the  Y-axis  a  line  of  sysaetry  of  the 
wave.  The  phase  velocity  C  is  defined  as  the  average  fluid  velocity  at  any  horizontal 
level  caapletely  within  the  fluid. 

We  introduce  a  potential  function  4(X,Y)  and  a  streaa  function  T(X,Y)<  Let  the 
stream  function  assume  the  values  zero  and  -Q  on  the  free  surface  and  on  the  bottom 
respectively.  The  undisturbed  fluid  depth  B  is  defined  by 

(2.1) 

we  take  the  origin  of  our  coordinate  systaa  on  the  undisturbed  level  of  the  free 
surface,  so  that  the  bottoa  is  given  by  Y  “  -H,  and  we  denote  the  equation  of  the  free  r 
surface  by  Y  -  ?(X). 

The  exact  nonlinear  equations  for  •  and  C  are 

0,  -H  <  Y  <  C(X>  , 


•«,  * 


■  0  on  Y  -  -H  , 
*X?X  "  *y  "  0  on  Y  "  C(X> 


(2.2) 

(2.3) 

(2.4) 


1  ,.2  .  ,2,  .  _  T 

2  <*X  +  V  +  9C  '  P 
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(1  ♦  V 


2\3/2“ib2  °n  Y-'(X> 


(2.5) 


■quation  (2.2)  follows  froa  conservation  of  Bass.  Equations  (2.3)  and  (2.4)  are  the 
conditions  that  the  bottoa  and  free  surface  are  streamlines,  and  equation  (2.5)  is  the 
dynamic  free  surface  condition.  The  Bernoulli  constant  b  on  the  right  hand  side  of  (2.5) 
is  to  be  found  as  part  of  the  solution. 

In  order  to  derive  asymptotic  solutions  it  is  convenient  to  introduce  the  following 
dimensionless  variables 


♦ 


(2.6) 


-3- 


In  (2.6)  L  la  a  length  scale  of  the  wave  in  the  X  direction  and  A  la  a  Measure  of  the 
amplitude.  Rewriting  (2. 2)-(2.5)  in  terae  of  the  dimensionless  variables  (2.6)  gives 


64  +4  -  0 

**  yy 


•1  <  y  <  an(x) 


4  -0  on  y  -  -t 

y 

*nx4x  ~  i  *y  “  0  y  *  m<*> 


(2.7) 

(2.8) 

(2.9) 


1  a*2  +  -  -  a2  +  n  -  8t _ SL_ 
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“  -r-  on  y  -  an(x)  . 


(2.10) 


In  (2.7)-(2. 10)  a  and  $  are  the  dimensionless  parameters 

A  H2 

«  -  f  “  2  '  (2.11) 

Relations  (2.7)-(2. 10)  are  the  classical  watar  wave  equations. 

Me  seek  solutions  for  periodic  water  waves  of  wavelength  X,  and  introduce  the 
dimensionless  wavelength 

t  -  X/L  .  (2.12) 

The  Froude  number  F  is  defined  by 


*“-“■*/  ♦<**  •  (2.13) 

✓5?  »  0  * 

Solitary  waves  are  the  limit  of  periodic  waves  as  i  ♦  «.  in  that  case  F  -  B. 

In  order  to  motivate  the  considerations  presented  in  the  next  sections,  it  is 
worthwhile  reviewing  briefly  some  of  the  classical  perturbation  solutions  of  the  system 
( 2.7)-( 2. 11 ).  Stokes  (1847,  1880)  derived  a  perturbation  solution  for  pure  gravity  waves 
by  assuming  a  small  and  0  of  order  one.  His  results  were  generalised  to  include  the 
effect  of  surface  tension  by  Harrison  (1909),  Wilton  (1915)  and  Mayfeh  (1970).  However, 
these  perturbation  calculationa  become  invalid  as  0  ♦  0  because  the  ratio  of  successive 
terns  is  then  unbounded. 

The  shallow  watar  equations  era  derived  by  assuming  0  small  and  a  of  order  one. 
These  equations  do  not  have  travelling  wave  solutions  because  the  dispersive  effects  are 
neglected  (Whithorn  (1974),  p.  4S7).  The  inclusion  of  dispersive  effects  into  the  shallow 


-4- 


water  theory  leads  to  the  Jtortawag  de  Vries  equation.  This  equation  can  be  derived  by 
assuming  0  small  and  a  of  order  0  (Keller  (1948),  vanden-Broeck  and  Shen  (1982)). 


Thus  we  let 


«  »  0  ■  * 


and  expand  n,  4  A1111  *  »■ 


n  -  «i„  ♦  «n,  +  e  n2  +  *•• 

♦  ♦  40  ♦  «♦,  ♦  «2#j  ♦ 

B  »  B0  ♦  CB1  *  •••  . 


(2.14) 


(2.15) 


»0  *  r0  "  1 
A1  at  «n 

fo  no 

2F,nl  "  3n.nl  ♦  (t  -  4>,'A" 
10  0  0  3  0 


Then  we  use  (2.14)  and  (2.15)  in  ( 2.7)-(2. 10),  and  find  that 

(2.16) 
(2.17) 

0  .  (2.18) 
The  primes  in  (2.17)  and  (2.18)  denote  derivatives  with  respect  to  x.  The  dispersive 
effect  arises  from  the  last  term  in  (2.18). 

Korteweg  and  de  Vries  (1895)  showed  that  periodic  solutions  of  (2.18)  can  be  found  in 
closed  form.  As  the  wavelength  tends  to  infinity  these  waves  tend  to  the  solitary  wave 
solution 

n„  -  a  sech2(|) 

(2.19) 

.  -  . 

When  T  <  1/3  these  are  elevation  waves  with  Fronds  number  greater  than  one  when 
x  >  1/3  they  are  depression  waves  with  Froude  number  less  than  one. 

Bquation  (2.19)  shows  that  the  slope  of  the  wave  profile  becomes  large  near  r  -  1/3 
and  the  solution  ceases  to  exist  altogether  when  x  -  1/3.  Thus  the  Ksrteweg  de  Vries 
equation  (2.18)  (denoted  in  the  remaining  part  of  the  paper  as  the  KdV  equation)  is  invalid 
in  the  neighborhood  of  x  -  1/3.  This  is  due  to  the  fact  that  the  dispersive  effects 
disappear  as  x  approaches  1/3. 


In  this  section  we  shell  derive  en  equation  analogous  to  the  MV  equation  which  is 
valid  in  a  neighborhood  of  r  -  1/3.  The  coefficient  of  the  dispersive  ten  nj'*  in 
(2.18)  vanishes  at  T  -  1/3  asking  a  travelling  wave  impossible.  To  obtain  a  balance 
between  the  dispersive  and  nonlinear  terns  near  t  »  1/3  we  take 


in  (2.7)-(2. 10).  Then  we  expand  n»  4#  B  and  t  as 

n  -  nQ  +  en,  +  e2n2  +  ••• 

4  —  ^2  +  4g  +  +  C^2  +  *^^3  +  *** 

(3.2 

1  2 

T  -  3  +  eTi  +  e  t2  +  ••• 

B  -  B0  ♦  eB,  +  eSj  +  ••• 

f  -  r0  +  ep,  +  c2r2  +  •••  . 

Prom  (3.1)  and  (2.7)-(2.8)  4  satisfies 

+  4  -  0 

“  ”  (3.3 

4  “0  on  y  -  -1  . 

y 

Then  using  (3.2)  in  (3.3)  and  equating  coefficients  of  powers  of  e  to  sero,  we  obtain 


♦0  “ 


♦  .--4(1  +  y)2f  +  g(x) 

1  a  1  i  ,  (3*4) 

42  -  77  d  +  y>4fuv)  --5(1  +  y)  V  +  h(x) 

43  -  -  |j-  (1  +  y)6f<vi>  +  77  <1  +  y)Vv  --5(1+  y>V  +  k(x)  . 

In  (3.4)  f,  g,  h  and  k  are  arbitrary  functions  of  x  and  prines  denote  differentiation 
by  x. 

We  use  (3.1),  (3.2)  and  (3.4)  in  (2.9)-( 2. 10),  expand  in  power  series  about  e  •  0 
0  1  2 

and  equate  coefficients  of  e  ,  e  and  e  .  This  gives  the  following  equations 

B0f  +  nQ  -  0  (3.5) 

B0n0  +  f"  "  0  (3.6) 


®«g'  +  n,  "  7  ♦  -5  »of'”  "  1 

Boni  ♦  95  ■  7  t(lv)  -  Vo 


®0*»*  ♦  *>2  ’  +  ^  9'"  "  4T  BOf<V>  -  V  -  \  (f’>2  *  i  ^  +  Vo  +  VaT  "  g,) 


Vi  +  h"  *  “Vo  +  6  9 


«  j.  1  _(*▼)  „  1  -(vi)  ...  | 

0  +  6  9  ■  V  ■  sT  f  ■  Vo  ■  Vi  • 


(3.10) 


For  (3.5)  to  have  a  nontrivial  solution  we  must  hava 


Bq  ■  *  • 


We  shall  consider  waves  wing  to  the  right  and  take 


B„  -  1  . 


(3.12) 


Then  (3.5)  implies 


(3.13) 


He  differentiate  (3.7),  add  the  result  to  (3.8)  and  use  (3. 12)— (3. 13).  Then  it 


follows  for  a  nontrivial  solution  that 


B1  “  0 


(3.14) 


Using  (3. 12)— ( 3. 14)  in  (3.7)  gives 


n.  -  -g*  +  4  f"'  . 


Next  we  use  (3. 12)~(3. 15)  in  (3.9)-(3. 10),  and  after  eliminating  n2  and  h  we  obtain  the 
following  equation  for  n. 


2B2ns  -  3vj  ♦  vr  -  h  C*  - 0  * 


(3.16) 


Finally  we  show  that  the  Bernoulli  constant  B  is  the  same  as  the  Proude  number  P 
to  the  order  considered  here.  From  (2.1),  (2.6)  and  (2.11)  we  have 


an 

F  “  «  /  ♦  <*y  • 


Then  using  (3.1),  (3.2)  and  (3.4)  in  (3.17)  we  find  that 

F„  -  B0  -  1,  Ft  -  B,  -  0,  Fj  -  B2  .  (3.18) 

h  similar  perturbation  calculation  using  (3.1)  can  be  performed  on  the  time  dependent 


water  wave  equation.  This  leads  to 


2(gHf  /2C  -  25 


x  -  i  «x  +  “2(t  -  T><xxx  -  h  CXXXXX 


Equation  (3.19)  reduces  to  (3.16)  for  travslling  wavs  solutions  5  «  5(X  +  ct). 

In  the  next  section  we  seek  periodic  solutions  of  (3.16)  by  a  Stokes  type  expansion. 


4.  Periodic  Waves  gwr  t  -  1/3 

He  seek  a  periodic  eolution  of  (3.16)  in  the  fora  of  an  expansion  in  powers  of  the 
wave  amplitude.  Therefore  we  write 

nfl  »  an01  +  *2n02  +  o(a3)  ,  (4.1) 

P2  -  Fjo  +  aF21  +  8^22  +  0(a3)  .  (4.2) 

Here  a  is  a  measure  of  the  wave  amplitude,  the  precise  meaning  of  it  will  become  clear 
later. 

Next  we  substitute  (4.1)  and  (4.2)  into  (3.16)  and  (3.18)  and  collect  all  terms  of 
like  powers  of  a.  Thus  we  obtain 

<  i.i 

0  . 


9V  vi*  +  T  n***  —  -  - 

"20n01  1n01  45  n01 


2r20n02  +  Tin02'  "  45  n02>  "  3n0in01 


2F  T) 1 
21  'oi 


The  solutions  of  (4.3)  and  (4.4)  are  given  by 


01 


'02 


cos  Kx 


-SKlIfiT^  -  8KF^ 


M  k5i-i 


45 


K  ]  'cos  2kx 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


20 


vl*k« 

2  90 


f2,  -  0 


(4.7) 

(4.8) 


Hera  X  -  2»/i  is  the  wavenumber.  Relations  (4.1)  and  (4.5)  define  a  as  the  amplitude 
of  the  fundamental  in  the  Fourier  series  expansion  of  nQ(x). 

The  classical  dispersion  relation  for  linear  water  waves  is  given  by  (see  Nhitham 
(1974),  p.  403) 

*2  -  -55  ♦  7  4^-)  .  (4.9) 

1/_ 

From  (2.11),  (2.12)  and  (3.1)  we  have  2IH/X  ■  e  '*  Ji/l,  Substituting  this  result  and 
(3.2)  into  (4.9),  and  expanding  for  e  small  we  find 

F2  -  1  +  2F20 

where  Fjq  is  given  by  (4.7).  Thus  the  solution  of  (3.16)  overlaps  the  classical  linear 
solution,  as  the  amplitude  tends  to  sero. 

The  perturbation  solution  (4.5)>(4.8)  is  invalid  when 


'20 


“  2t  K2  +  —  r* 
1*  45  * 


(4.10) 
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4 


because  tha  coefficient  of  cos  2Kx  in  {4*6)  is  then  unbounded*  When  (4*10)  is  satisfied/ 
the  solution  of  (4.3)  is  given  by 

■  cos  Kx  +  EjCos  2Kx  (4.11) 


'20 


TiK  JeI  2T  t2  +  —  r4 

2  90  *1*  45  ^  ' 


(4.12) 


Here  b2  is  a  constant  to  be  found  as  part  of  the  solution.  Substituting  (4.11)  into 


(4.4)  we  obtain 


1  ( v) 

2v  t*  +  t  n" •  -  ~  n 
**20  0  2  ln02  45  n02 


(»21K  -  -  n2)sin  Kx 


+  (4P21KE2  -  ^  K)sin  2Kx  -  ~  KEj  sin  3Kx  -  3K(B2)2Sin  4Kx  .  (4.13) 

The  solutions  of  (4.13)  are  periodic  and  bounded  if  and  only  if  the  coefficients  of 


sin  Kx  and  sin  2Kx  in  the  right  hand  side  of  (4.13)  vanish.  Thus  we  have 


(4.14) 


P21  *  ±  \  —  • 

21  4  /2 

Substituting  these  results  into  (4.1)  and  (4.2)  we  obtain 

nQ  -  a  cos  Kx  ±  s2~  ^2  COS  2Kx  +  0(a2)  , 

2 


(4.15) 


(4.16) 


V 


7 .  w 

2  2 


±  3a2”5/2  +  0(a2)  . 


(4.17) 

Relations  (4.16)  and  (4.17)  show  that  two  solutions  exist  when  (4.10)  is  satisfied. 
Similar  properties  were  found  by  Wilton  (1915),  Pearson  and  rife  (1961),  Schwarts  and 
Vanden-Broeck  (1979)  and  Chen  and  Saffman  (1980)  for  waves  in  water  of  infinite  depth  and 
by  Nayfeh  (1970)  for  waves  of  small  amplitude  in  water  of  moderate  depth.  It  is 
interesting  to  note  that  the  solutions  (4.13),  (4.14)  given  by  Nayfeh  (1970)  converge  to 
(4.16)  and  (4.17)  as  the  depth  tends  to  xero.  Thus  Nayfeh's  solution  and  the  present  long 
wave  calculation  overlap. 

A  non-uniqueness  analogous  to  the  one  described  by  (4.16)  and  (4.17)  will  occur  in 
general  when  waves  of  wavenumber  K  and  nK  travel  with  the  same  phase  velocity,  i.e. 


when 


(4.18) 


v.  V 

2  90 


n2t  K2  4  4 

1  n  K 

2  90 


Hera  n  is  an  arbitrary  integer  greater  than  one.  The  solution  of  (4.3)  is  then 

n0  “  cos  Kx  ♦  Bn  cos  nKx  . 

For  n  •  2,  (4.18)  and  (4.19)  reduce  to  (4.11)  and  (4.12). 

Relation  (4.18)  can  be  rewritten  as 


(4.19) 


45T. 


1  ♦  n 


2  ' 


(4.20) 


It  follows  f row  (3.2).  (3.14)  and  (4.20)  that  this  non-uniqueness  can  only  occur  for 
T  <  1/3. 

Numerical  values  of  Fj  for  a  -  0.01  and  K  ■  1  obtained  by  integrating  (3.16) 
numerically  are  shown  in  Figure  1.  The  broken  line  corresponds  to  the  solution  (4.7), 
(4.8)  and  the  two  crosses  to  the  solutions  (4.17).  These  results  indicate  that  the  two 
solutions  (4.17)  are  members  of  two  different  families  of  solutions. 

Other  numerical  solutions  of  ( 3.16)  will  be  discussed  in  Section  7. 
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Values  of  F2  as  a  function  of  r  for  a  «  0.01  and  K  -  1 .  The  solid  curves 
correspond  to  the  numerical  integration  of  (3.16),  the  broken  line  to  the  solution 
(4.7),  (4.8)  and  the  two  crosses  to  the  solutions  (4.17). 


Reformulation  as  an  Integro-dlf ferentlal 


itlon 


It  la  convenient  to  reformulate  the  problem  aa  an  integro- differential  equation  by 
considering  the  complex  velocity  u  -  iv.  Hare  u  and  v  are  the  horisontal  and  vertical 
components  of  the  velocity  respectively.  The  variables  are  made  dimensionless  by  using 
H  as  the  unit  length  and  C  aa  the  unit  velocity.  He  choose  the  complex  potential 

f  -  ♦  +  i*  (5.1) 

as  the  independent  variable. 

In  order  to  satisfy  the  boundary  condition  (2.8)  on  the  bottom  f  -  -1,  we  reflect 
the  flow  in  the  boundary  ♦  -  -1 .  Thus  we  seek  x  +  ty  and  u  -  iv  as  analytic  functions 
of  f  in  the  strip  -2  <  ♦  <  0. 

Next  we  define  the  dimensionless  wavelength  t  “  A/H  and  introduce  the  change  of 
variables 

f  -  ♦  +  i*  -  -i  log  p  (5.2) 

where 

16 

P  ”  re  .  (5.3) 

Relations  (5.2)  and  (5.3)  map  the  bottom  f  -  -1,  the  free  surface  t  *•  0  and  its  image 

2» 

Sf  9 

♦  ”  -2,  respectively,  onto  the  circles  r  “  r0  -  e  ,  r  -  1  and  r  «  rjj. 

The  values  of  the  real  and  imaginary  parts  of  the  function  G(p)  «  u  -  iv  -  1  on  the 
free  surface  r  -  1  and  its  image  r  -  r^  are  related  by  the  identities 

u<9)  -  1  Real{G(ei0)}  »  1  +  Real{G(r*ei8) }  ,  (5.4) 

v(4)  -  -Im{G(e18)}  -  Im(G(r*ei8)}  (5.5) 

where 


The  functions  u(9)  and  v(9 )  in  (5.4)  and  (5.5)  are  the  horisontal  and  vertical 
components  of  the  velocity  on  the  free  surface  9  ■  0. 

In  order  to  find  a  relation  between  u  and  v  we  apply  Cauchy's  theorem  to  the 
function  G(p)  in  the  annulus  r*  <  I p I  <  1.  Using  the  relations  (5.4)-(5.6)  and 
exploiting  the  bilateral  symmetry  of  the  wave  about  9-0  we  find  after  some  algebra 


5  1/2 

u(0)  -  1  -  -  f  /  (u(s>  -  1]da 

*  n 


t/2 


jf  v(s) [cotg  7 


irJ  / 

l  0  i 


Ft/2  [u(a)  - 


+ 


2 

t 


Ft/2  [u(a)  - 


(a  -  6)  ♦  cotg  7  (a  +  6)]da 

1J[r2  -  coa  (B  -  0)]  ♦  v(a)ain  &  (a  -  0) 
1  ♦  t*  -  2r2eoa  (a  -  0) 

1]  [r2  -  coa  -2|  (a  +  6)1  +  v(a)ain  (a  ♦  0) 
1  r*  -  2r2coa  ■—  (a  +  0) 


da 


da  . 


The  third  integral  in  (5.7)  ia  of  Cauchy  principal  value  fora. 
The  aurface  condition  (2.10)  can  now  be  rewritten  aa 


1  2  2  1  2  2 

\  F^u^r  +  \  F^v^r  +  / 


v(a) 


0  [u(a)]2  +  [v(a)l 2 


ds 


-  T  -  1  f2[u(o,] 

Uu(0)]2  +  (»(♦))  2>  '2 


2  Tv1 ( 0 ) 
u(0) 


In  the  remaining  part  of  the  paper  we  ahall  chooae  coordinates  x,y  with  the 
at  a  creat  or  a  trough  of  the  wave.  The  ahape  of  the  free  aurface  la  then  defined 
paraaeterically  by  the  relations 

i(+)  -  u(a){[u(a)l  2  +  (v(s)]2}“1da  , 

0 

y(0)  -/  v(e){(u(a)l  2  +  (v(a)]2}  'da  . 

0 

Finally  we  inpose  the  periodicity  condition 

x(  t/2)  -  t/2  . 

He  shall  measure  the  amplitude  of  the  wave  by  the  paraaeter 

Uq  -  u(0)  . 


(5.7) 


(5.0) 

origin 

(5.9) 

(5.10) 

(5.11) 

(5.12) 
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For  given  values  of  t,  Ug  and  (,  (5.7)-(5.12)  define  a  systee  of  integro- 

differential  equations  for  u(+)<  v(|),  *(♦),  y  ( ♦)  and  r. 

The  equations  for  solitary  waves  are  obtained  by  taking  the  Unit  i  ♦  •  in  (5.7) 
This  leads  after  iaae  algebra  to 


«ce>-i v<,,lr~^s  ♦  rrrlds  *  z  I  *  -  6->-y('>  il  „ 

0  evsve  w0  (■  -  8)2  ♦  4 


+  If  i£-t  »W>  ±  zJl  a.  . 

*  0  (a  ♦  8)J  ♦  4 


C5.1 


In  the  next  section  we  describe  a  mis r leal  achene  to  solve  these  equations 


(a)  Par iodic  wavea 


To  aolva  tha  ayataa  (5-7)— (5.12)  we  introduce  tha  M  mesh  points 

. *•  <«•»> 

We  also  define  tha  corrasponding  quantities 

Uj  -  u(4z),  1  -  1.....N  .  (6.2) 

V1  "  1  “  ,«***»"  •  (6.3) 

It  follows  frost  symmetry  that  v1  «  v^  -  0,  so  only  H  -  2  of  tha  Vj  ara  unknown.  Xn 
addition  fro*  (5.12)  wa  have  Uj  -  Uq  so  only  N  -  1  of  the  u^  ara  unknown.  We  shall 
also  use  tha  N  -  1  midpoints  y  given  by 

♦x*V2  "  2(*x*W'  1-1 . .  -  1  ‘  (6‘6> 

wa  evaluate  uI+  y2  -  u(*t+^),  ^\2~  v<*1+ 1/2>-  4+ V2  “  "’‘♦i+V2>'  “*d 

Vj+  y^  ■  v* (4I+  y  )  by  four  points  interpolation  and  difference  formulas. 

wa  now  diacretise  (5.7)  by  applying  the  trapesoidal  rule  to  tha  integrals  on  the 
right-hand  aide  with  s  *  4^,  X  -  and  6  -  $I+  y^  ,  I  »  1,  H  -  1.  The  symmetry  of 

the  quadrature  formula  and  of  tha  discretisation  enables  us  to  evaluate  the  Cauchy 
principal  value  integral  as  if  it  were  an  ordinary  integral.  In  this  way  we  obtain 
N  -  1  algebraic  equations. 

Next  we  substitute  into  (5.8)  the  expressions  for  1/2 .  yi+V2*  “i+  V2  *nd  vi+ V2 
at  tha  points  4  ..  ,  I  “  2, ...,N  -  1.  The  integral  in  (5.8)  is  evaluated  by  the 

l+  '/j 

trapesoidal  rule  with  the  mesh  points  a  ■  The  derivative  v'(0)  in  (5.8)  is 

approximated  by  a  four  points  difference  formula.  Thus  we  obtain  another  N  -  2  algebraic 
equations. 

For  given  values  of  t,  Ug  and  i  we  have  therefore  2N  -  3  algebraic  equations  for 
the  2N  -  2  unknowns  u^,  Vj.  and  F.  The  last  equation  is  obtained  by  discretising 
(5.11). 


The  2N  -  2  equations  are  solved  by  Newton's  method.  After  a  solution  converges  for 

V 

given  values  of  T,  uQ  and  t,  the  surface  profile  x(  ♦) ,  y(>)  is  obtained  by  applying 
the  trapesoidal  rule  to  (5.9)  and  (5.10). 


For  Mch  calculation  presented  in  Section  7  the  number  of  mesh  points  H  was 
progressively  increased  up  to  a  value  for  which  the  results  were  independent  of  M  within 
graphical  accuracy.  Most  of  the  computations  were  performed  with  M  ■  60.  A  few  runs  were 
performed  with  N  -  100  as  a  check  on  the  accuracy  of  the  computations . 

In  the  remaining  part  of  the  paper  we  shall  refer  to  the  above  mmerical  scheme  for 
periodic  waves  as  numerical  schasm  I. 

(b)  Solitary  waves 

To  solve  the  aystam  (5.8)— (5.13)  we  introduce  the  N  mesh  points 

♦j  -  S(I  -  1),  I  -  1,...,K  .  (6.5) 

Here  E  is  the  interval  of  discretisation. 

The  quantities  Uj,  Vj,  ♦I+1/2'  <»i+  1/2  *nd  vJ+  y2  are  defined  as  in  the  previous 
subsection.  It  follows  from  symmetry  that  vt  -  0.  In  addition  u,  -  Ug  so  that  only 
2N  -  2  of  the  Uj  and  Vj  are  unknown. 

The  discretisation  of  (5.13)  is  entirely  analogous  to  the  procedure  used  to  diecretise 
(5.7).  In  this  way  we  obtain  H  -  1  algebraic  equations.  The  truncation  error  due  to 
approximating  the  Infinite  integrals  by  integrals  over  a  finite  range,  was  found  to  be 
negligible  for  ME  sufficiently  large.  As  shown  in  the  previous  subsection,  the 
discretisation  of  (5.8)  leads  to  another  N  -  2  algebraic  aquations. 

Thus,  for  given  values  of  T,  uQ  and  i  we  have  21-3  algebraic  equations  for 
the  21  -  1  unknowns  Uj,  Vj.  and  P.  The  last  two  equations  are  obtained  by  imposing 
-  1  and  vH  -  0. 

We  shall  refer  to  this  numerical  scheme  as  numerical  scheme  II. 


7.  Discussion  of  the  Results 

(a)  Depression  solitary  waves 

numerical  schesn  II  was  us  ad  to  coaiputa  daprasslon  solitary  wavas  for  given  valuas  of 
t  and  Ug.  In  ths  first  calculation,  tha  Nawton  i tar at Iona  wara  started  with  tha 
asymptotic  solution  (2.19)  as  the  initial  guess,  for  Ug  close  to  one,  the  iterations 
converged  rapidly.  Once  a  solution  was  obtained  it  was  used  as  the  initial  guess  for  the 
next  calculation  with  a  slightly  dlffsrent  value  of  T  or  Ug.  The  curve  a  in  Figure  2 
shows  a  typical  profile  for  t  *  0.7  and  Ug  ■  3.0. 

The  asymptotic  solution  (2.19)  can  be  rewritten  in  terms  of  the  variables  used  in  the 
numerical  scheme  as 


y  ■  h[sech2  [■ 


3A 


4(1  -  3t) 


),/2~  1 
1  *-1] 


(7.1) 


F2  -  1  ♦  A  .  (7.2) 

The  curve  b  in  Figure  2  corresponds  to  the  profile  (7.1)  in  which  the  amplitude  A  is 
equal  to  the  amplitude  of  the  numerical  solution.  For  1  <  Ug  <  1.03  the  numerical 
results  and  the  asymptotic  formula  (7.1)  were  found  to  be  indistinguishable  to  graphical 
accuracy. 

In  Table  1  we  compare  numerical  values  of  F  with  the  approximation  (7.2)  for  various 
values  of  uQ  and  T  -  0.4,  and  0.7. 

As  Ug  increases  for  a  given  value  of  t>  j,  the  wave  profile  becomes  steeper  and 
the  distance  between  the  trough  and  the  bottom  decreases.  For  t  >  this  distance  tends 
to  zero  as  uQ  ♦  <*.  The  corresponding  Froude  number  tends  to  zero  and  the  profile 
approaches  a  static  limiting  configuration  in  which  gravity  is  balanced  by  surface 

tension.  Then  (2.10)  reduces  to  a  differential  equation  for  the  free  surface: 

2  "3/2 

n  -  70^(1  nR)  -  1  .  (7.3) 

The  boundary  conditions  for  (7.3)  are 

0(0)  -  0  (7.4) 

0(»)  -  1  .  (7.5) 
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TABLE  1 


“0 

T  -  0, 

nunarical 

4 

KdV 

T  •  0 

numerical 

.7 

KdV 

1.0 

1.0 

1.0 

1.0 

1.0 

1.03 

0.986 

0.986 

0.986 

0.985 

1.5 

0.844 

0.850 

0.832 

0.831 

2.0 

0.756 

0.764 

0.734 

0.734 

4.0 

0.578 

0.590 

0.541 

0.542 

10.0 

0.396 

0.407 

0.35S 

0.355 

20.0 

0.291 

0.301 

0.255 

0.253 

50.0 

0.190 

0.195 

0.163 

0.152 

100.0 

0.136 

0.137 

0.116 

0.081 

Table  It  Values  of  the  Froude  number  for  depression  solitary 
waves  for  t  «  0.4  and  T  »  0.7,  and  various 
values  of  uQ.  Numerical  values  were  computed  by 
scheme  II,  and  the  KdV  values  found  from  (7.2)  with 
A  taken  equal  to  the  amplitude  of  the  numerical 
solution.  * 
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Multiplying  (7.3)  by  nx  and  Integrating  with  reapect  to  x  yields 


-2  5  "1/2  i 

—  +  t(i  ♦  nx)  -  n  +  T  -  -  . 

The  value  of  the  constant  of  integration  in  (7.6)  was  evaluated  by  using  (7.5). 
Integrating  (7.6)  gives  a  formula  for  the  shape  of  the  free  surface/  namely 


(7.6) 


x  -  /  [t2(t>  ♦  t 

0 


1 

2 


2 


(7.7) 


The  curve  c  in  Figure  2  represents  the  profile  (7.7)  for  T  «  0.7. 

Me  denote  by  tan  0  the  slope  of  the  profile  (7.7)  at  the  trough  x  •  0.  Then  (7.4) 
and  (7.6)  yield 

cos  0  «  1  -  .  (7.8) 

Relation  (7.8)  implies  that  the  present  limiting  configuration  is  only  possible  for 

x  >  1/2. 

For  T  <  1/2,  the  numerical  computation*  indicate  that  the  wave  ultimately  reaches  a 
critical  configuration  with  a  trapped  bubble  at  the  trough.  This  critical  configuration  is 
shown  in  Figure  3  for  t  -  0.4.  similar  limiting  configurations  were  obtained  previously 
by  Crapper  (1957),  Schwartz  and  Vanden-Broeck  (1979)  and  Vanden-Broeck  and  Keller  (1980). 
Haves  for  larger  values  of  uQ  could  be  obtained  by  allowiong  the  pressure  in  the  trapped 
bubble  to  be  different  from  the  atmospheric  pressure  (Vanden-Broeck  and  Keller  (1980)). 

Numerical  scheme  II  was  used  to  compute  depression  solitary  waves  for  t  <  1/3.  in 
Figure  4  we  present  solutions  for  Uq  -  1.03  and  various  values  of  t.  As  t  decreases 
the  profiles  develop  a  large  number  of  inflexion  points.  He  were  unable  to  compute 
solutions  for  uQ  »  1.03  and  T  <  0.21  because  too  many  mesh  points  were  required. 

In  Figure  5  we  compare  the  numerical  solution  of  the  exact  nonlinear  aquations  with 
the  profile  obtained  by  numerically  integrating  (3.16).  He  found  that  the  two  solutions 
become  identical  within  graphical  accuracy  in  the  limit  as  t  ♦  1/3  and  uQ  ♦  1  with  the 
ratio  (Ug  -  1 ) ( x  -  1/3 )“ 2  constant.  This  constitutes  an  important  check  on  the 
consistency  of  our  results. 
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igure  3a  -  Computed  free  surface  profile  of  a  depression  solitary  wave  for  t  »  0.4 
and  Uq  »  227.  The  free  surface  has  just  one  point  of  contact  with 
itself  and  encloses  a  small  bubble  at  the  trough. 


Figure  3b  -  The  bubble  of  Figure  3a  expanded  by  a  factor 
of  12.5.  The  vertical  scale  is  the  same 
as  the  horizontal  scale. 
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Computed  free  surface  profile  of  a  depression  solitary  wave 


Figure  4  c  -  Same  as  Figure  4a  with  uQ  «  1.03,  t  “  0.228 
F  -  0.954. 


(b)  elevation  solitary  waves 

In  preliminary  computations  we  attempted  to  calculate  solitary  waves  for  T  <  1/3  by 


using  numerical  scheme  II.  The  iterations  converged  rapidly.  However,  the  profiles  were 
not  independent  of  M  and  S  and  they  did  not  appear  to  approach  a  limiting  profile  as 
MX  ♦  •  and  x  ♦  0.  Furthermore  at  large  x  the  profiles  oscillated  with  finite  amplitude 
about  the  undisturbed  level  instead  of  approaching  a  uniform  stream. 

On  the  other  hand  we  were  able  to  compute  periodic  waves  of  large  wavelength  t  for 
T  <  1/3,  by  using  numerical  scheme  I.  A  large  number  of  dimples  are  present  on  these 
profiles  (see  Figure  6).  We  found  many  different  families  of  periodic  waves.  This  non- 
uniqueness  agrees  with  the  perturbation  calculation  of  section  4.  waves  in  different 
families  are  characterised  by  the  number  of  dimples  in  a  wavelength.  In  Figure  7  we  plot 
the  Froude  number  versus  wavelength  for  3  such  families,  when  T  -  0.24.  Curve  a 
corresponds  to  14  dimples,  curve  b  to  15  dimples,  and  curve  c  to  16  dimples.  Graphs  of 
representive  members  of  each  of  these  families  are  presented  in  Figure  6. 

As  can  be  seen  from  Figure  7  each  family  only  exists  for  a  limited  range  of 
wavelengths.  In  particular  there  is  a  maximum  value  of  the  wavelength  beyond  which 
solutions  in  a  given  family  cease  to  exist.  Therefore,  a  solitary  wave  cannot  be  obtained 
as  a  limit  of  solutions  belonging  to  a  single  family. 

In  conclusion  all  our  attempts  to  find  elevation  solitary  with  surface  tension  by 
either  numerical  scheme  I  or  II  failed.  Mathematically  the  question  of  the  existence  of 
these  elevation  waves  remain  open. 


-27- 


4 


# 


Coaputed  tree  surface  profile 


QUf6  6C  -  Same  as  in  Figure  6a  with  1  -  36.6 
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